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We present a measurement scheme capable of achieving the quantum Hmit of parameter estima- 
tion using an adaptive strategy that minimizes the parameter's variance at each step. The adaptive 
rule we propose makes the scheme robust against errors, in particular imperfect readouts, a critical 
requirement to extend adaptive schemes from quantum optics to solid-state sensors. Thanks to 
recent advances in single-shot readout capabilities for electronic spins in the solid state (such as Ni- 
trogen Vacancy centers in diamond), this scheme can be as well applied to estimate the polarization 
of a spin bath coupled to the sensor spin. In turns, the measurement process decreases the entropy 
of the spin bath resulting in longer coherence times of the sensor spin. 



A common strategy for estimating an unknown param- 
eter associated with a field is to prepare a probe and let 
it interact with the parameter-dependent field. From the 
probe dynamics it is possible to derive an estimator of 
the parameter. The process is repeated many times to 
reduce the estimation uncertainty. A more efficient pro- 
cedure takes advantage of the partial knowledge acquired 
in each successive measurement to change the probe-field 
interaction in order to optimize the uncertainty reduc- 
tion at each step. This adaptive Bayesian estimation 
strategy has been proposed to improve the sensitivity 
of parameter estimation in quantum metrology [1]. It 
has been shown that adaptive estimation can achieve the 
Heisenberg or quantum metrology limit (QML) without 
the need for entangled states [2-6]. Here we introduce a 
novel adaptive scheme that attains the QML, as mani- 
fested by various statistical metrics of the estimated pa- 
rameters. In addition, the proposed scheme can be made 
robust against errors so that the QML is achieved e.g. 
even for imperfect readouts, a critical requirement to ex- 
tend adaptive schemes from quantum optics to solid-state 
sensors. We further present an application of the adap- 
tive scheme to the measurement of a quantum parame- 
ter: given single-shot readout capabilities for electronic 
spins in the solid-state [7-9], the scheme could be used 
to create a narrowed state of a surrounding spin bath, 
thus increasing the sensor coherence. In this context, the 
QML scaling translates into a shorter time for the nar- 
rowing process, an important feature when dealing with 
a finite bath relaxation time. 

Consider a two-level system {|0), |1)} interacting with 
an external field characterized by the parameter 6, "H = 
baz- A typical situation is a sensor spin-i interacting 
with a magnetic field. The parameter can be estimated 
by a Ramsey experiment (Fig. 1), where the probability 
of the system to be in the \m) state (to = {0, 1}) at the 
end of the experiment is given by 



VAm\b) = ^[1 



(-l)™e-"/^=cos(6T-Hi?)] (1) 



where ^ is the phase difference between the excitation 
and readout pulses and we introduced a decay with a 
constant T2 during the interrogation time r. If we have 
a prior knowledge of the parameter -described by an a 
priori probability distribution (p.d.f.) P^'^'{b)~ the mea- 
surement updates our knowledge, as reflected by the a 
posteriori probability: P{b\m) oc P^'^' {b)V^{m\b) . 

More generally, after each measurement we can update 
the probability for the phase ip = br, so that after n such 
measurements with outcomes to„, we have a p.d.f. 

P(")(<^|TO„) cx p("-i)((^|m„_i)7'4TO„|<^) (2) 

Thanks to the periodicity of the probability P{(p), we can 



expand it in Fourier series [4], P^^^ip) 
that we can rewrite Eq. (2) as 
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The proportionality factor is set by imposing that Pq — 
-^ as required for a normalized p.d.f. We can further 
generalize this expression when the system is let evolve 
for an integer multiple tn of the time r, thus obtaining a 
general update rule for the p.d.f.: 



(n) 1 (n-1) 

Pk « 2Pfc 



-UT/T2 , 



„i(m„7r+i5„)„("-l) i p-i(m„7r+)?„)„("-l) 
^ Pk-tn '^^ Pk+t„ 



(3) 



An adaptive strategy will then seek to choose at each 
step the optimal tn and i?„ that lead to the most efficient 
series of N measurements for a desired final uncertainty. 
In order to design an adaptive strategy, we need to 
define a metric for the uncertainty (and accuracy) of the 
estimate. The Fourier transform of the p.d.f. can be 
used to calculate the moments of the distribution as well 
as other metrics and estimator. From the formula for the 
moments, {(p"-) = J_ P{(p)(p°'dip = J^kP'^ I- e'^'^'^ ip" dip 
we can calculate the variance, 



* pcappell@mit.edu 



(ip^) - {iff = — po 
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where the average is {ip) — —'2inJ2kM k~Pk- 

The variance is often not the best estimate of the un- 
certainty for a periodic variable [3]. A better metric is 
the Holevo variance [10], 



y^ = (27r|(e^^)|)-2-l = (27r|p.i| 



1, 



(4) 



where we used the fact that (e*''') — p.i. We further 
notice that while the absolute value of p_i gives the phase 
estimate uncertainty, its argument provides an unbiased 
estimate of (p. More generally, estimates are given by 
ipest = arg((e**'^))/i = arg(p_i)/t, giving a new meaning 
to the Fourier coefficients of the p.d.f. 

The goal of the estimation procedure is then to make 
\p-i\ as large as possible. Assume for simplicity ip — 
br = and neglect any relaxation. Then the probability 
of the outcome m^ = is 7'i;(0|0) = i(l - cos^). We 
assume that we do not have any a priori knowledge on 
the phase, so that P('''((p) = 1/2tt. We fix the number 
of measurements, iV, each having an interrogation time 
T„ = tnT = 2^-"r [4, 11, 12]. A potential strategy would 
be to maximize |pi\ | at each step n. However, under the 

assumptions made, p_\ 
where it is 



until the last step, n = N, 
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Writing p_2 = qe^^, we have 



47r|pl/| = \/l + 47r2(j2 + 47rgcos(x + 2i9 



N) 



This is maximized for z9jv 
by maximizing q 



IpL^-^^ 



-x/2 = 5arg(pL^ ^^) and 
A similar argument holds 
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for the maximization of \p_2 \'- one has to set i!} 

I arg(p_4 ) and maximize |p_4 |. By recursion we 
have that at each step we want to maximize 



l/tl = 
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We have thus found a good adaptive rule, which fixes 
t„ = 2^-" and ^„ = i arg {p^Z'!)- 

With this rule we obtain the standard quantum limit 
(SQL) for the phase sensitivity, as we now show. Using 

the optimal phase, the Fourier coefficients p_/ 
each step 



are at 
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Then, for a total number of measurements N, the Holevo 
variance is Vh ^ {1 ~ 2-(^+i))-2 _ i « 2"^. The total 
interrogation time is T = t(2^+^ — 1) yielding 



y«m - ^. 



At 



(5) 



We can improve the sensitivity scaling and reach the 
QML by a simple modification of this adaptive scheme. 
Instead of performing just one measurement of duration 
tn at each n*'' step, we perform two, updating the p.d.f. 
according to the outcomes. For (/? = the update rule at 
each step is now 



(") 1 



6Pfc 



^Pfc-f„ +^Pk+t„ + Pk-2t„ + Pk+2t„ 



with the normalization factor 



(n) 

Restricting the formula above to the terms p_^ gives 
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By recursion this yields 



\P. 
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from which we obtain a Holevo variance that follows the 
QML, Vh ~ 3-2^^^, or in terms of the total interrogation 
time 



Vh 



A8Tt^{T + At) 
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(r-2T)2(T + 6T)2 r2 



(7) 



The classical and quantum scaling of the adaptive 
scheme with one or two measurements per step is con- 
firmed by the p.d.f. obtained in the two cases (Fig. 1). 
For one measurement, the final p.d.f Fourier coefficient 
are \pk\ = ^^(1 — 2~(^+^)|fc|) and the probability is well 
approximated by a sine function. 



F(^)(^) 



■)N+1 



27r 



-sinc(2^^+V) 



which gives a variance a ~ 2 ^' ^. For two measurements 
per step, instead, the p.d.f. is well approximated by a 

73 n-N 

2 ■ ^ 



Gaussian (see Appendix) with a width a 




-0.03 -0.02 -0.01 0.01 0.02 0.03 
cp (radians) 



-1000 -500 500 

k(radians-i) 



FIG. 1. P.d.f (left) and its Fourier transform (right) after 
an 8-step adaptive measurement, with 1 (red, dashed) and 
2 measurements per step (black). In the inset, Ramsey se- 
quence. 



We now consider possible sources of non-ideal behav- 
ior. The first generalization is to phases (p 7^ 0. In 
this case, while the SQL is still achieved with the one- 
measurement scheme, two measurements per step do not 
always reach the QML. Indeed, at each step there is a 
probability V{l — V) that the two measurements will give 
different results; if this happens at the n}^ step, we ob- 

(n) 

tain p_( = 0, thus failing to properly update the p.d.f. 
While the probability of failure is low, a solution could 
be to perform three measurements and update the p.d.f. 
only based on the majority vote. 

We can further consider the cases where the signal de- 
cays due to relaxation or there is an imperfect readout. 
Then the probability (1) becomes 

Vi>{m\b) = i[l-c(-l)"e-^/^^cos(6T-|-z9)], 

with c the readout fidelity. Considering the effects of only 
this (constant) term, the update rule Eq. 6 becomes 
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FIG. 2. Holevo variance vs. total time T = Mr(2'^+^ - 1), 
with T — 1. (•) with 1 measurement per step (M=l) Vh 
follows the SQL. (o) 2 measurements per step (M=2) 
achieve the QML. (*) with c=0.95, the QML scaling is 
lost, but can be preserved for longer with M=4 (a) and 
restored (n) by setting M=n + 1, even for lower c, e.g. 
c=0.85 (•) 



p-L 



(n-l) 



1 



+ c^p. 



(n-l) 



We can calculate a recursion relationship in the limit of 
good measurement, e = (1 — c) « 0, to obtain 



\p-i 



1 

2^ 



1 



-(l + iVe)2- 



-2N 



which yields an Holevo variance Vh ~ 3(1 -I- eN)2^'^^ 
that does not follow anymore the QML scaling, except 
for eiV '--^ 1. A similar, more complex result is expected if 
relaxation effects are taken into account (see Appendix) . 
A strategy to overcome this limitation is to repeat the 
measurement at each step more than two times (Fig. 2). 
Specifically, setting the number of measurements M=n -I- 
1 (if allowed by relaxation constraints) restores the QML 
scaling. 

The proposed adaptive method promises to achieve 
Heisenberg-limited estimation of a classical phase with- 
out the need of fragile entangled states and thus it could 
improve the sensitivity e.g. of recently proposed mag- 
netic sensors [13]. It can be as well used to measure 
a quantum variable, such as a phase resulting from the 
coupling of the sensor to a larger quantum system or 
bath. In turns, the measurement can be used to lower 
the entropy of the bath (usually a thermal equilibrium 
mixture) yielding an increase in the coherence time of 
the sensor [11, 14, 15]. The QML scaling of this adap- 
tive method translates into a faster narrowing of the bath 
dispersion, which would improve similar schemes in solid- 
state systems [16, 17], where the bath itself might present 
fluctuations. 

Specifically, we consider the coupling of a sensor spin to 
a spin bath. This situation is encountered in many physi- 
cal systems, such as quantum dots [18, 19] or phosphorus 
donors in silicon [8, 20]. Here we analyze as an example 



the system comprising a Nitrogen- Vacancy (NV) center 
electronic spin coupled to the bath of nuclear ^^C spins 
in the diamond lattice [21, 22]. Recent advances in the 
measurement capabilities [7] offer single-shot read-out of 
the NV state, thus enabling adaptive schemes. 

In a large magnetic field along the NV axis, the hy- 
perfine interaction between the electronic spin and the 
nuclear spins is truncated to its secular part, % = 
•^z X]fc ^klz,k — SzAz (where S denotes the electronic 
spin, Ik the nuclear spins). During a Ramsey sequence 
on resonance with the nis = 0,1 energy levels of the elec- 
tronic spins, the coupled system evolves as 

m)) = [sin {A,t) ]1) + cos {A,t) ]0)] ]V)c, (8) 

where \^)c is the initial state of the nuclear spin bath. 
The measurement scheme (Ramsey followed by NV read- 
out) is a quantum non-demolition measurement [23-25] 
for the nuclear spins, since their observable does not 
evolve - as long as the secular approximation holds. 
The adaptive process is then equivalent to determin- 
ing the state-dependent (quantized) phase ip — (Azt). 
The uncertainty on the nuclear bath state, pc = 
^aPa lV'a)(V'alc' ^^ reflected in the p.d.f of the phase 
(with an injective relation if the operator Az has non- 
degenerate eigenvalues). Thus updating the phase p.d.f. 
will update the density operator describing the state of 
the nuclear bath. After each readout of outcome to, the 
system is in the state 



p(") ex \m){m\p^"-^^\m){m\ 

= \m){m\Y,a'Pi}{m\ipa)Pa \lpa){'4'a\ 



(9) 



withPtf(TOl(/?„) = l(TO,Val[sin(A,i)|l)+COs(^,i)|0)]|V'a)| 

Note that in this expression the probability update rule 
is equivalent to Eq. 2 and thus the adaptive procedure 
ensures that the final state has lower entropy than the 
initial one. 
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FIG. 3. Simulated NV center spectrum from a Ramsey se- 
quence before (dotted) and after (solid line) an 8-step adap- 
tive measurement. The adaptive scheme narrows the bath 
distribution, here corresponding to a 50 x improvement of the 
coherence time. We considered ~ 2600 closest nuclear spins, 
randomly positioned in the diamond lattice, in a IT magnetic 
field and initially in the maximally mixed-state. We simulated 
full dipolar couplings among the nuclear spins and between 
the NV and the bath spins (approximating the hyperfine in- 
teraction) to validate the secular approximation. The narrow 
spectrum reveals an average nuclear field of (B^) « 700nT 
(this was chosen at random among possible nuclear spin state- 
dependent fields; similar results were obtained for different 
fields and spin positions in the lattice). In the inset, Holevo 
variance (circles) as a function of the total time, and signal 
decay (dotted line) due to intra-bath couplings during the 
adaptive scheme interrogation time (TRamscy)- 



A difference between measuring a classical field and a 
quantum operator is that in the latter case the resulting 
phase is quantized, thus it has a discrete p.d.f.. An ex- 
treme case is when all the couplings to the Nc nuclear 
spins are equal, Ak = a, V/c. Then the eigenvalues are 
na/2, with \n\ < Nc integer, each with a degeneracy 
d{n) = (^ /2+n)- While the adaptive scheme needs to 
be modified (e.g. by considering a discrete Fourier trans- 
form) , we note that since all the eigenvalues are an integer 
multiple of the smallest, non-zero one (a for Nc even, a/2 
for Nc odd), we only need M steps, with 2^ > ^ [with 
minimum interrogation time r = 27r/(a2*^)], to achieve 
a perfect measurement of the degenerate phase tp [11]. 
In the more common scenario where Ak varies with the 
nuclear spin position (and Nc is large enough) the eigen- 
values give rise to an almost continuous phase [15], thus 
it is possible to directly use the adaptive scheme derived 
above. 

As an example of the method, we consider one NV 
center surrounded by a bath of nuclear spins (^'^C with 
1.1% natural abundance). At low temperature and for 



NV with low strain, it is possible to perform single-shot 
readout of the electronic spin state with high fidelity in 
tens of /is [7] . Optical illumination usually enhances the 
electronic- induced nuclear relaxation [26] , due to the non- 
secular part of the hyperfine interaction. This effect is 
however quenched in a high magnetic field {B > IT) and 
the relaxation time is much longer than the measurement 
time (Ti > 3ms [27]), sign of a good QND measurement. 

We simulated the Ramsey sequence and adaptive mea- 
surement with a bath of ~ 2600 spins around the spin 
sensor in a large magnetic field. We considered the full 
anisotropic hyperfine interaction between the NV and the 
^^C spins and we took into account intra-bath couplings 
with a disjoint cluster approximation [28, 29]. Even for 
the longest evolution time of the Ramsey sequence re- 
quired by the adaptive scheme, the fidelity F of the sig- 
nal with the ideal Ramsey oscillation (in the absence of 
couplings) is maintained. After an 8-step adaptive mea- 
surement, the nuclear spin bath is in a narrowed state. 
We note that in general the adaptive scheme does not po- 
larize the spin bath (indeed a final low polarization state 
is more probable). However, the bath purity is increased, 
which is enough to ensure longer coherence times for the 
sensor spins, since it corresponds to a reduced variance 
of the phase and hence of the sensor spin dephasing. In 
Fig. 3 we compare the NV center spectrum for an evo- 
lution under a maximally mixed nuclear spin bath and 
under the narrowed spin bath. The figure shows a re- 
markable improvement of the NV coherence time. 

In conclusion, we described an adaptive measurement 
scheme that has the potential to achieve the quantum 
metrology limit for a classical parameter estimation. We 
analyzed how imperfections in the measurement scheme 
affect the sensitivity and proposed strategies to overcome 
these limitations. This result could for example improve 
the sensitivity of spin-based magnetometers, without re- 
curring to entangled states. In addition, we applied the 
scheme to the measurement of a quantum parameter, 
such as arising from the coupling of the sensor to a large 
spin bath. We showed that the adaptive scheme can be 
used to prepare the spin bath in a narrowed state: as the 
number of possible configurations for the spin bath is re- 
duced, the coherence time of the sensor is increased. The 
scheme could then be a promising strategy to increase the 
coherence time of qubits, without the need of dynamical 
decoupling schemes that have large overheads and inter- 
fere with some magnetometry and quantum information 
tasks. 
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Appendix A: Probability distribution for the adaptive scheme with two measurements per step 

In the main text we considered the adaptive scheme where at each step two measurements (with the same reading 
time) are carried out and presented an approximate formula for the p.d.f that is obtained after N steps. Here we 
present details of the derivation for the case where ip = 0. 
In the Fourier space, the p.d.f. coefficients can be calculated from the recursive relation to be 



27r 



(2~+^-l-|fc|)(2"+^-|fc|)(2~+^ + l-|fc|)-4(2~+^-l-|fc|)(2»+^-|fc|)(2»+^ + l-|fc|) , . „JV+1 „ o 
Pl^ ^ ^ "" (2" + ^-l-|fc|)(2"+^-|fc|)(2~+^ + l-|fc|) ^ ^ 2JV+1 _ 2 

fc < 2^+1 - 2 
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For large N we can simplify the expressions as: 
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In turn, these expressions are well approximated by a Gaussian (although the original function has longer tails). 
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Appendix B: Influence of noise on the adaptive scheme 



In the main text we discussed the effects on the adap- 
tive scheme of imperfect readout of the sensor state. A 
similar effect is expected as well if the signal decays due 
to decoherence during the Ramsey interrogation time, 
as the difference in the probability of getting a different 
result (m — 0, 1) given a different phase is reduced by 
a factor e""^^"^^ : 



ri)il\b) ~ 7'^(0|6) = le-^/^' cos(5r 



^). 



The time-dependence of such an imperfection, though, 
complicates the recursive relationship (Eq. 3 of the 
main text), thus we analyze the effect of decoherence 
numerically. We find that the decay effectively sets a 
maximum number of measurements (see figure), since 
the interrogation time cannot exceed T2. A small im- 
provement is achieved by increasing the number of 
measurements per step, but the QML is not recovered 
at longer times. 

Another source of imperfection would derive from 
variations of the phase during the total estimation 
time. The rate of this variations sets an upper limit 
to the number of steps in the adaptive scheme. 
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FIG. 4. (color online) Holevo variance vs. total time in the 
presence of signal decay. Black lines with circles: M — 2 mea- 
surements per step. Red lines with squares: M=n -I- 1 measure- 
ments at the n"* step. The relaxation constant T2 was taken to 
be T2/T = 5 X 10"* (dotted line), 10"^ (solid lines), 2.5 x 10"^ 
(dash-dotted lines) and 5 x 10""^ (dashed lines). The total time 
T was measured in units of the dimensionless time r. 
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